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ABSTRACT 


An attempt lo made to validate a method that uoes radiometric 
surface temperatureo and a boundary layer model to eotimate ourface 
energy budgets and characteristics. The surface temperatures are 
area-overaged values so the results from the Radiometric Data/Model 
System (RDMS) are also area-averages. Another remote sensing method, 
sodar. Is used to obtain heat flux estimates on the some scale as 
the RDMS values for ground truth measurements. A simultaneous 
collection of radiometric surface temperatures from a hand-held 
radiometer and sodar data was made on seven days between mid-July 
and mid-October 1980. The comparison of the ROMS and sodar heat 
fluxes proved disappointing. Free convection conditions, required to 
producr sodar-derived heat fluxes, were Inhibited by a terrain- 
induced low level Inversion. Only three out of seven cases produced 

meaningful sodar heat fluxes. Of those three cases, one had good 

-2 

agreement and the other two had sodar heat fluxes 15 to 45 W m 
lower than the RDMS values. Since the RDMS method is relatively 
untested, it was impossible to conclusively determine Its validity 
from the results. There was evidence that the true heat flux was 
not underestimated by the RDMS, so It could be concluded that the 
Bowen ratios over well-watered vegetation were likely to be quite 


small. 
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1.0 INTEODUCTION 

1.1 DEVELOPMENT OP THE SATELLITE DATA/MODEL OVSTEM (SDMS) 

Carloon and Boland (1973) dovolopod a ono^dlmonoional boundary 
layor modol capable of prodieElnQ Eho ourfaeo EomporaEuro and Eho 
ourfaeo onoroy budgoE. They uood Eho modol to oKamlno Eho rolaElvo 
ImporEaneo of Eho variouo Eorrain paromoEoro govornlng Eho rooponoo 
of Eho ground Eo heaElng. From a oorloo of oonolEivlEy EooEo, 

Ewo paramoEoro, Eho ground moloEuro availablllEy and Eho ground 
conduccancQ (Ehormal InorEia), woro found Eo domlnaEo Eho rooponoo 
of Eho ground Eo oolar radlaElon. Thuo» Ehey f@lE Eha£ if Eho 
valuoo of moloEuro avoflablllEy (M) and Ehormal InorEia (P) for a 
ourfaeo could bo obEainod, iniEializaEion of Eho oEhor variabloo 
in Eho modol uoing oEandard mo Eooro logical daEa oourcoo would produco 
ground EomporaEuroo accuraEo Eo ± 2G. 

Thoy noEed EhaE blade body ourfaeo EomporaEuroo clearly ohov;ed 
Eho rooponoo of a given Eype of ourfaeo Eo hoaEing. Thio rooulE 
suggeoEod a way of doEermining valuoo of M and P for a ourfaeo. 
Carlson and Boland proposed Eo combine an invoroion of Eho modol 
ouEpuE 1 ‘fiEh a pair of observed radiomoEric surface EemporaEuroo . 
Surface EemperaEure variaEions such ao Ehooo measured by a oaEolliEe 
would be maEched to Ehe variaEions calculaEed by Ehe modol. From 
EhaE maEching, unique values of M and P for a surface would be 
defined. Given esElmaEes of Ehooe Ewo domlnanE paramoEero, the 
model could then calculate Ehe heat flux and evaporation for Eho area 
observed by Ehe satellite. Thus, Carlson and Boland felt that Eho 
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ouboEFaeo ehoraceorioeieo and ourfaeo onorcy budooe for an area 
could bo ooEiraoeod from oaEolllEo oboorvoElono coQbincd wish Eho 
ouEpuE from ’i ono^dimonoional boundary layor modol. 

AuouoElno (1978) accoloraEod Eho linkino of oaEolliEo daEa tijiEh 
Eho modol by dovolopino a Eochnlquo Eo produce ourfaeo EomporoEuro 
mapo for an area uolng radiomoErle oaEolliEo daEa. He producod day 
and nighE Eomporaturo mapo of Loo Angoloo. UnforEunaEoly* Eho 
oaEolliEo overpaoo Elmoo woro juoE afEer ounrloo and ounooE, and 
Eho mapo produced did noE ohow Eho raonlmum diurnal EomporaEuro 
variation. The matching tochnlquo propoood by Carloon and Boland 
roqulroo that ourfaeo tomporaturoo bo moaourod near Eho Elmoo of 
maximum and minimum heating. Thuo» no attempt wao made In AuguoEino'o 
work to Implement their procedure. 

It wao not until 26 April 1978 that a oaEolliEo capable of 
producing effective black body tomperature moaouremonto of the earth 
txooe to the times of maximum and minimum tomporaturoo v;ao launched. 

On that day, a small applications -'xploror oatolllte wao launched 
by the National Aeronautics and Space Administration (NASA) as part 
of the Heat Capacity Mapping Misolon (HDMM) The operational 
characteristics of the HCMM satellite are discussed by Barnes and 
Price (1980). By chat time, the method for linking the onc-dlmcnslonal 
boundary layer model with satelllte*>sensed surface temperatures had 
been perfected- 

When HCF I data became available, the satellite data/model system 
was used to produce maps of M, P, heat flux, and evaporation for 
St. Louis and Los Angeles. Carlson et al. (1981) detailed the results 


oS ehooo eoooot Frea tho ©opo @£ M and P» choy ©ado e@neiuQi©no 
eoneoffiino eho ©oehanioao behind eho urban eoaperoturo anoaaly. Tho 
mapo of hooe flun and ovaporoeion poreroyod fho opaeial variation 
of eho pareieionino of oolar radiaelon cauood by difforonooo in land 
ucto ond vocotativo eovor, 

Addieional uork wao done uoino HCtCl data by Kocin (1979) and 
DiCrioeofaro (1980). CCocin onaainod eho paecorno of M, P, hoat fltn: 
snd Gvoporaeion over a vocotatod watorohod In Miooouri. Ho wao obio 
to relate opociite rogienol footuroo (forooto, eroplando, ote.) and 
rainfall to the analyoio of moioturo. availability. Changoo in tho 
diotrlbution '>f moioturo avoilabllity wore found to bo ee^dominatod 
by rainfall amounto ai'd changoo in tho ourfaec'o vogotativo canopy. 
DiCrlotofaro produced mape of tho ground charaetoriotieo ond energy 
budget for both urban and rural aroao. Variationo in moioturo 
avoilabllity and thermal inertia wore ohoxm to bo rooponoible for 
the temperature d.'fferencoo between urban and rural roglono. Hio 
rooultd tor the rural areao indicated thot fho variability in heat 
fiujc and evaporation wao caubod by variationo in land uoo. 

1.2 STATEMENT OF THE PROBLEM 

Before employing the satellite data/model oyotom (SDMS) for 
fn. ther investigations, the validity of ito products needed to bo 
established. Although the previouoly mentioned invootigatoro have 
demonstrated the ability of the oyotem to produce reasonable patterns 
of the surface energy budget and substrate characteristics, tho valuoo 
of these parameters needed confirmation by independent measurements. 
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UnforeunatGly, independent moaouremento of the oamo nature and 

type as the model products are not easily obtained. The estimates 

of all the parameters produced by the SDMS customarily use radiometric 

2 

surface temperature measurements averaged over approximately 1 lun . 

It follows that the values of moisture avollability, thermal inertia , 
heat flux, and evaporation are also area°averaged . A large array of 
point measurements is usually required to determine area<=>averaged 
ground truth values for these parameters. However, the financial 
and physical resources for such a measurement program can be enormous. 
Thus, here we have used other methods in an attempt to verify the 
SDMS results. 

One such method is aircraft measurements. DlCristofaro (1980) 
attemt-ted to verify some of his results utilizing this procedure. 

He obtained aircraft measurements of surface heat and evaporative 
fluxes from the Sulfur Transport and Transformation in the Environment 
(STATE) project (Schlermelr et al., 1979). A comparison between 
results from the SDMS and the STATE data was attempted. Unfortunately, 
the results were inconclusive. The agreement was reasonable for the 
evaporative fluxes, but the aircraft measured values of heat flu.: 
were sevcial times larger than the satellite derived values. It was 
noted, however, that the aircraft measured large heat fluxes in a 
forested area which displayed low daytime temperatures on the satellite 
pictures. In fact, the aircraft-measured heat fluxes for this area 
were only slightly lower than the values found in downtown St. Louis. 

Another remote sensing method, sodar, can produce values of heat 
flux averaged over the same scale as the satellite measurements. The 
sodar measures the scattering of sound in a volume of air by temperature 
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fluctuations. The scattering can be related to the temperature 

2 

structure constant C^. l-Jhen free convection condltioVf.*j e:tiet in 

the atmosphere, the surface heat flux obeys a 4/3 power lav; betv;een 
2 

and height (Wyngaard et al., 1971). In this situation, the largo 

2 

values are concentrated in updrafts generated by the heating of 

the terrain and advected by the mean wind. Neff (1973) found regions 
2 

of large that were localized at a horizontal scale of about 300m. 

2 

Heat flux estimates derived from values during free convection 
conditions were therefore averaged over a section of the terrain 
upwind of the sodar antenna. The exact area of the average depended 
on the mean wind speed. 


1.3 PURPOSE OF THE THESIS 

A need existed to verify the values of the SDMS by Independent 
measurement. One technique that can measure heat flux on the same 
scale as the SDMS is sodar. This research compared sodar measure- 
ments of heat flux with those from a version of the SDMS that used 
surface temperature measurements from a hand-held radiometer instead 
of a satellite. Radiometric temperature and sodar data were collected 
simultaneously on seven days in the summer and autumn of 1980 to 

provide the raw data for the two methods. The sodar data was 

2 

processed by producing ten minute averages of profiles. If free 
convection conditions were present in the atmosphere, an estimate of 
the heat flux for that period could then be made from the profile. 

The heat flux estimates of the modified SDMS and the sodar were 
analyzed over the diurnal cycle. An effort was made to determine the 
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validity of the modified SDMS values by comparing them with the sodar 


results. 
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2.0 METHODS OF ESTIMATING THE HEAT FLUX 

2.1 THE SATELLITE DATA/MODFi. SYSTEM (SDMS) METHOD 

The onQ-dlmenslonal boundary layer model used and Che procedure 
of linking the model with the satellite data in the SDMS have been 
well documented (see Carlson and Boland (1978), Carlson et al. (1981), 
Dodd (1979), Kocin (1979), and DiCristofaro (1980)). Thus, only a 
brief summary of the SDMS is presented here. 

The purpose of the SDMS is to produce maps of substrati 
characteristics and surface fluKes for an area using surface tempera- 
ture measurements near the maximum and minimum of the heating cycle. 
The maps are created by matching the obseirved surface temperatures 
to surface temperatures calculated for the observation time by the 
model. The basic assumption which allows the transformation between 
surface temperatures and substrate parameters to be made is that a 
given model maximum/minimum temperature pair is computed from a 
unique pair of M and P values and vice-versa. Thus, a measured 
temperature pair can be combined with the model output to infer 
specific values of M and P. Carlsoa and Boland (1978) found that M 
and P are the dominant factors in determining the partitioning of 
incoming solar radiation into sensible and latent heat flux at the 
surface. Accordingly, solutions yielding a unique M and P also 
correspond to a unique solution of the surface energy budget. The 
model output at the time of the temperature measurements for a wide 
range of M and P values is represented by a set of regression 
equations Involving the maximum/minimum temperature pairs as 


predictors. In this manner, oubotrato characteristics and the surface 
energy budget are obtained when an observed majclmum/mlnimuiD temperature 
pair is Inserted Into the regression equations. Manlmum/minlmum 
surface temperature maps from an aircraft or satellite can be trane° 
formed using the regression equations Into maps of the t^i;o surface 
parameters and the surface fluxes. 

2.2 SODAR HEAT FIUX THEORY 

A monostatic sodar measures acoustic energy backscattered by 

the atmosphere. This energy can be related to a temperature structure 
2 

parameter which, in turn, can be used to estimate heat flux. 

The radar equation for backscattered energy can be modified to 
apply to a monostatic sodar. Following Unden^ood (1981), a form of 
the monostatic sodar equation Is 

P(r) = Pq “I E a(ir) (2.1) 


where 


P(r) = received acoustic power, 

P^ = transmitted acoustic power, 

A == cross-sectional area of the antenna, 
r 

r = range of scatterer from antenna, 
a = path averaged molecular attenuation coefficient, 
= length of the scattering pulse in space, 

E ° turbulence induced excess attenuation, 
a(7T) = backscattering cross section. 


9 

The backocacterlng croQO section o is the fraction of the 
incident power backscattered per unit solid angle. With proper 
calibration of the sodar system and estimates of the attenuation, 
Equation 2.1 can be used to obtain values of o. Those values are the 
basis for estimating heat flux with the sodar. However, the calcula- 
tion is done indirectly through another variable, the temperature 

2 

structure parameter C^, which is related to both a and heat flux. 

2 

is defined as 

° (T(x) - T(x4-r)]/r^^^ (2.2) 

2 

where represents the variance of the temperature difference 

between two points a distance r apart and r is assumed to be in the 

inertial subrange. Using Tatarskl (1971) it can be shown that a 
2 

is related to in this fashion: 

a = 0.0039 k^'^^ . (2.3) 


Wyngao'.d et al. (1971) derived the free convection limit 
-2 _ 3 ,T.2/3 ,"o,4/3 ,-4/3 

where 

k ° the von Karman constant, 

T = the average air temperature, 
g = gravity, 
p = density, 


(2.4) 
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C ° tho specific heat at constant pressure, 

P 

H ° the heat flux, 
o 

Coulter and Wesely (1980) added a humidity correction factor Co 
equation 2.4 v;here 

° Hi + (0.07/S)^]/(l + (0.06/0)^)}^^^ 

SO 

X (1 + 0.07/0^)’^ . (2.3) 

Thus, when free convection conditions exist, sodar, by measuring a, 
can be used to estimate the heat flux. 
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3.0 DATA COLLECTION 

3.1 DESCRIPTION OF FIELD EXPERIMENT 

The SDMS requires measuremonco of the ourfoce tomporaturo near 
the imuclmum and minimum of the heating cycle to produce cotimates 
of heat flux. The SDMS v;as designed to utilize temperature data 
from polar orbiting satellites that pass over a site ot the proper 
points in the heating cycle. The surface temperatures from these 
satellites represent an average value for an area of one pixel, 
which is about one half to one square kilometer for the highest 
resolution infrared satellite radiometers. SDMS values of heat 
flux were Intended originally to be calculated from satellite data. 
Unfortunately, no satellite capable of producing surface temperature 
measurements at the optimum times was operational during the period 
of the experiment. Thus, it v;as necessary to simulate satellite 
measurements by conducting a small field experiment to collect the 
temperature data required to calculate the SDMS heat flux values. 

The necessary surface temperature values near the maximum and minimum 
of the heating cycle vjere recorded using a hand-held radiometer in 
an observation technique, to be described below, which was designed 
to determine an area-averaged temperature value. This modified 
version of the SDMS was called the Radiometric Data/Model System (RDMS) . 
The acoustic signal returns of the sodar were collected at the same 
time and stored on magnetic tape to later be converted into heat 
fluxes. 

The locale for the field experiment was the Rock Springs 
Agricultural Research Center of The Pennsylvania State University. 
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This facility in oituatod oovoral miloo from tho Ponn State cotnpuo 
(oQG Figure 3.1) and waa chooon becauao it io tho oito of tho 
Meteorology Departmont'o eodar antenna. The oroo ourrounding the 
Qodar conaiatQ mainly of heterogoneoua cropland, in v;hlch largo 
fielda of alfalfa, corn, and ooyboano pluo omall toot ploto of 
alfalfa, corn, and oata aro prooont within 300m of tho antonna. A 
picture of the Rock Springs site is prooontod in Figure 3.2. 

The duration of the experiment was July to October 1980. Since 
the boundary layer model employed in the ROMS procedure can simulate 
only clear sky and non*-advectlve conditions, the temperature and 
sodar data were collected only when those conditions were present. 

In all, seven days met the criteria for data collection as demonstrated 
by the strong diurnal variation in the air temperature traces recorded 
at Penn State for each case (see Figures 3. 3-3. 9). A summary of the 
dates and weather conditions of each data collection day is presented 
in Table 3.1. 


3.2 SURFACE TEMPERATURE DATA COLLECTION 

The radiometer used in this research was a Barnes PRT-5 on loan 
from the Environmental Protection Agency (EPA). The radiometer was 
calibrated at Penn State by using a calibration box where the 
radiometer views a 'black box' cavity with a water bath of variable 
knotm temperature. The temperature reading of the radiometer Is 
compared with the water temperature. It was discovered through this 
process that the radiometer consistently gave readings 1.5 to 2.0C 
too warm. Thus, all temperature readings from the field were corrected 
before final use by subtracting 1.7C from the raw value. 
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Figure 3.2 A section of the Rock Springs Agricultural Research 
Center of The Pennsylvania State University where 
the sodar and surface temperature data for this 
experiment were collected. 
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Ths idea of uolng a hand^^hold radiomotsr (aloo called an Infrared 
thermometer) to make surface temperature measurements of vegetation 
has been exploited previously on numerous occasions. Fuchs and 
Tanner (1966) found that, given the emlssivlty of the surface and the 
reflected atmospheric radiation, they could measure vegetative canopy 
temperatures with a radiometer to ± O.IC. 

Fuchs et al. (1967) measured crop surface temperatures using a 
radiometer. They examined the dependence of the temperature on viewing 
angle and solar azimuth angle (see Figure 3.10). 

The Incident angle had a minor effect on the measured temperature 
v;ith variations usually less than 1C. Lotror temperatures were 
recorded for Incidenr angles near normal and greater than 60 degrees. 

At the small viewing angles, the radiometer looks deeply into the 
vegetation and sees more of the shaded (cooler) portions of the 
vegetation and ground surface. At large incident angles, the radiometer 
sees plant tips plus the horizon and therefore can produce a lower 
reading. 

The solar azimuth angles had an effect In the case of Fuchs et 
al. (1967) only when the radiometer was viewing row crops. In this 
case, when the radiometer was facing the sun and the crop rows formed 
an angle greater than 10 degrees with the sun, temperatures were 
consistently 1.3C lower than when the radiometer viewed the same 
crops while pointed away from the sun. This difference was confirmed 
by leaf temperature measurements made with thermocouples. The 
radiometer views primarily shaded leaves at a solar azimuth angle of 
0 degrees, versus the primarily sunlit leaves when the solar azimuth 
angle of the radiometer is 180 degrees. Fuchs et al. also found a 
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Figure 3.10 Scheniaclc of the incident angle a» and the solar 
azimuth angle (j> of the radiometer sensing beam. 
The field of view of wide and narrow beam Instru 
ments is given by and 62 * respectively. 

(Source: Fuchs et al,, (1967)). 
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± .3C variation in time when the viewing angle wao kept eonotant. 

They felt thlo variation wao cauood by air turbulence changing the 
ventilation of the crop leaveo, A 30 degree field of view radiometer 
produced Icoo variation in temperature meaouromento than a 7 degree 
field of view instrument oinee it averaged the temperature over a 
larger area* 

The observation technique for the radiometric surface temperatures 
in the Rock Springs data collection v;as designed to enable one 
temperature to be determined for the entire site. A hand-held 
radiometer v;lth a narrow field of view can only scan a small area. 

Thus, to produce an average value for the site, the observation 
technique had to satisfy two objectives: 

1. Maximize the area that the radiometer viewed in one scan. 

2. Provide a large number of temperature readings to obtain 
temperature values of all the crops present at the site. 

The first step in maximizing the area viewed by the radiometer 
was to use a wide angle field of view lens. Further, the radiometer 
was operated at a height of about 5m above the surface from a site 
atop a tower mounted on a mobile van (see Figure 3.11). Finally, 
temperature measurements were made at three or more locations by 
driving the van/ tower platform around the site. 

Twelve temperature observations were made at each site. Because 
of the surface temperature variations with incident angle and solar 
azimuth angle noted by Fuchs et al. (1967), the observations were 
made at three incident angles at each of four solar azimuth angles 
(see Figure 3.12). The temperature values were recorded both manually 
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Figure 3.11 


Van/ tower platform that surface temperature measurements 
were made from using a hand~held radiometer. 


0 “ 



Figure 3.12a The four solar azimuth angles at 
which surface temperature data 
were measured. 
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and aueoiaaeically by tnoano of a otrlp chare recorder , Addielonally, 
eho type of vocotatlon viewed by eho radleaoeer for each oboervaelon 
wao noted, ao v;ore ehe local weather condltlono (ooo Table 3.1). 

Unlike Pucho and Tanner (1966), in thio field ouperimont the 
vosotativo emiooivlty woo aoouBod to bo equal to one, althouQh plant 
canopioo typically have an emiooivity of about .97. Therefore, the 
amount of atmoopheric dovniitJollino longwave radiation rofloetod by 
the vegetation wao not calculated. The high degree of accuracy that 
Pucho and Tanner achieved with their radiometric temperature 
meaauremento wao not required for thio o::periment. A detailed 
diocuooion of the error produced by aoouming a vegetative emiooivity 
of one io given in Appendix A. The reoulto ohow that the tempera tureo 
measured in thio experiment are accurate to ± 0.3C. 

Thus, the result of the radiometer part of the data collection 
was 36 surface temperatures measured over a variety of vegetation 
types. Temperature observations were made twice on a data collection 
day; first at approximately 1200 EST and again Just after 0000 EST 
the following night. 


3.3 SODAR DATA COLLECTION 

The sodar used in the Rock Springs data collection was a 
modified Aerovlronment Model 300. The sodar was calibrated on June 11, 
1980, using the method described by Underwood (1978). No post- 
experiment calibration was performed. 

Th^ received acoustic energy of the Aerovlronment is usually 
qualitatively recorded on a facsimile display. At Penn State, however. 
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eho oyoEGiQ hoo boon modiflod £© aliea eho roeoivod aicnol to bo both 
quQneiCQtlvoly dioplayod in faloo color Imacory ond loQOod on 
masnoeic eapo. 

FIbupo 3.13 lo a block dlacrasn of the oodar oyotoa uood in ehio 
onporimont. Tho Dota Gonoral Novo 2/10 ainicoiopueor triOBOfO a burot 
I'' 1600 hs acoustic onorgy froa tho Aorovlronotont evory olght ooeonda. 
Aftor a pauoe of .3 mo» the oodar bocoisoo a rocolvor and convorto 
the backocatterod ocouotic onorgy into a voltago. Tho voltage oignal 
io tranoformed in a coupler for tranomlooion via phono lino from 
the oite back to tho computer laboratory. Thoro a demodulator 
convorto tho oignal back to a voltage that io oomplcd at 30 hz by an 
analog to digital converter. Tho Novo 2/10 procoosoo tho digltiaod 
signal, displaying it in faloo color on a video terminal and otoring 
it on magnetic tape. 

This entire proceoo is automated. On tho morning of a data 
collection day, the data logging program woo entered into tho 
minicomputer and the phone link t-;lth Rock Springs established. 

Then the sodar data t*;as logged until the evening of that day at which 
point the data collection program wau terminated by the user. 
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4.0 DATA ANALYSIS 

4.1 INTRODUCTION 

After Che last data collection day on October 19» 1980. the 
process of converting the raw sodar and temperature data Into heat 
flmces began. The temperatures were analyzed to produce a single 
day/nlghc pair of surface temperatures for each case. The temperature 
pair v;as then used In the RDMS to determine the heat flux for that 
case. The sodar data was plotted and examined for periods when free 
convection conditions vtere present. During those periods, heat 
fluxes were estimated with the sodar. 

4.2 SURFACE TEMPERATURE DATA 

The temperature data collection process yielded 36 day and night 
temperature values for each case. Figures 4. 2-4. 8 display the raw 
temperature data for each case and Che symbols are explained In 
Figure 4.1. Note the large amount of variation present. The average 
difference for the seven cases between the high and low daytime values 
is IOC. Highlighted In Figures 4. 2-4.8 are the three main causes of 
this variation: 

1. The viewing angle of the radiometer. 

2. The solar azimuth angle of the radiometer. 

3. The type of crop viewed by the radiometer. 

A further discussion of each of these factors follows. 

Fuchs et al. (1967) showed that when using a 30 degree field of 
view lens, extreme incident angles (> 60 degrees) tend to produce 
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PAY TEMPERATURES 


Solar Azimuth Angle 


0 " 

90° 

180° 

270° 

70° oza 

C3Z3 


C2Z3 

Incident 50° 

QQ3 



Angle 350 


EM 2 a 



NIGHT TEMPERATURES 


Incident 

Angle SO^DS 

35°czza 


Vertical lines Inolde a boK 
Indicate two or moro measure^ 
ments of that value were 
obtained. 


CROP ABBREVIATIONS 


C ° Corn G ° Grass 

A = Alfalfa S ° Soybeans 


SHADED BARS 


The shaded bars represent the best estimate 
of two standard deviations around the area- 
averaged temperature determined from the data. 


Figure 4.1 Explanation of the symbols used in the 
radiometric surface temperature charts 
(Figures 4, 2-4. 8 ). 
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cooler readings. An examination of the 70 degret incident angle 
temperatures (marked with a single o*" Figures 

4.2“4.8) shows that these temperatures v/ere among the coolest for 
every crop type. The average 70 degree incident angle temperature 
was 1“2C cooler than the average 50 degree temperature and 1=4C 
cooler than the 33 degree temperature. It is clear, since the 
radiometer was equipped v;lth a 20 degree field of view lens, that at 
an incident angle of 70 degrees the temperature reading from the 
radiometer vlev;ed a portion of the horizon. 

The temperature difference in row crops viewed at different 
solar azimuth angles found by Fuchs et al. vras also present in the 
Rock Springs temperature data. Corn was the only row crop at the 
Rock Springs site, and the average corn temperature at a solar 
a-imuth angle of 0 degrees (marked by and '///' in Figures 

4. 2-4. 8) is about 2C cooler than the average value measured at other 
azimuth angles. 

The spread in the temperature data was also Increased by the 
variety of crops at the Rock Springs site. Whether a particular 
crop is green and growing or brown and senecent causes a different 
reaction to the incident solar radiation. Radiometric measurements 
made by Blad and Rosenburg (1976) and Heilman and Kanemasu (1976) 
showed that leaf temperatures can vary between crops under the same 
micrometeorological conditions due to differences in evapotranspira- 
tion between the canopies. The grass at Rock Springs was characterized 
as "dry" in observations from all seven cases because of its brown 
color. Note in Figures 4. 2-4. 8 that the grass temperatures were 
consistently among the highest measured. When the soybeans were 
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well watered and green, as in Cases 2 and 3 (Figures 4.3 and 4.4), 
their temperatures were some of the lov/est found. In Cases 4 and 3 in 
September, the soybean field v;as characterized as "discolored (yellow) 
and only 25% green." This effect was caused by lack of watering 
during a very dry period. In the earlier cases the soybean tempera” 
tures were in the lower and middle portions of the temperature distribu- 
tion. In Case 6 (Figure 4.7), the soybean temperatures were as high 
as the dry grass values. At this time (late September), the entire 
crop had dried out. Finally, in Case 7 the soybeans were harvested 
and only a dried stubble remained. The radiometer was then essentially 
viewing dry ground and the soybean temperatures were the highest 
values recorded in Case 7. 

The large variability in the temperat* res caused by the three 
aforementioned factors made it difficult to come up with an area- 
averaged value. A subjective analysis procedure was used that 
involved examining the temperatures as displayed in Figures 4. 2-4. 8 
and usir three guidelines: 

1. The average temperature at an incident angle of 70 degrees 
was 1 to 4C cooler than the temperatures from the 50 and 
35 degree incident angles, 

2. It is the crop temperatures upwind of the sodar antenna 
which determine the sodar measured heat fluxes. 

3. The more area a crop covered, the greater its contribution 
to the heat flux. 

The temperature value and uncertainty for Case 1 (the shaded 
bar in Figure 4,2) provide an example of this process. Note that all 
but 2 of the 11 temperatures omitted at the low end of the error bar 
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are from a 70 degree Incident angle. The temperatures at the other 
extreme were eliminated for two reasons: 

1, Although the soyba?ns covered a large area, they v;ere 
located downwind of the sodar antenna. 

2. Grass covered a very small portion of the site, and therefore 
would not ha."' had a l^.r^e effect. 

This subjective analysis procedure yielded a temperature value 
with an uncertainty that was felt to be equivalent to two standard 
deviations. The shaded bars in Figures 4.3-A.8 are the estimated 
temperatures and uncertainties for Cases 2-7. The subjectively 
determined temperatures were within .5 to .15C of the average 
temperature of the cr^^ps upwind of the sodar antenna (alfalfa, corn, 
and grass) in all seven cases. Additionally, the average 50 and 35 
degree incident angle temperatures of every case were within 1C of 
the subjective value or within the calculated uncertainty. 

The nighttime temperatures are displayed in Figures 4. 2-4. 8 in 
somewhat the same manner as the daytime values, except there was no 
need to specify the solar azimuth angle. The average spread in the 
nighttime temperatures is only about 4C. The view factor problem of 
the 70 degree values is still present, but the vegetation has had 
enough time to reach a fairly uniform temperature. The nighttime 
temperature value and error for each case was determined through the 
same subjective analysis method as the daytime temperatures. However, 
because of the small spread of the night temperatures, the nighttime 
uncertainty was smaller than the daytime version. Note that no 
nighttime temperature data was obtained for Case 5. Thus, the value 
determined for Case 4 (the previous night) was used in the ROMS. 
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At this point, the process of calculating heat fluxes with the 
ROMS began. First, the boundary layer model was initialized for each 
case with data obtained from standard meteorological sources such 
as facsimile maps. Table 4.1 lists values of some of the Important 
input variables. 

The model vras also initialized with a sounding that is a smoothed 
version of the morning (12Z) observation. During July and August 1980, 
radiosondes were being launched every six hours here at Penn State as 
part of the Northeast Regional Oxideat Study (NEROS) . These local 
sounding data were used in the first three cases. For Cases 4 through 
7, the soundings from Pittsburgh (about 20( km away) were used. Since 
the data collection days were chosen for their lack of advection, the 
Pittsburgh data were considered representative of the local conditions. 

Next the model was put through a 'productiot run' where a wide 
range of surface temperatures and fluxes for the same initial 
conditions were created by varying the values of moisture availability 
(M) and thermal Inertia (P). The output from a 'production run' of 
the model was used to create a set of regression equations. These 
regression equations inverted the model output and thus were able 
to predict values of M and P for every case from the observed day 
and night temperatures. The predicted M's and P's were then used 
to calculate the RDMS heat fluxes. The model was run again for 
each case using the same set of initial conditions, but M and P 
were fixed at their predicted values. The heat flux curve from 
this model run was taken as the RDMS heat flux estimate for a given 
case. Appendix C contains an example of the process for determining 


a heat flux estimate. 


Table 4,1 Model input parameters for each case 
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In order to project the uncertainty in the temperature values 
through the ROMS, day and night temperatures from either end of 
their respective error bars were used to create tv/o additional day/ 
night temperature pairs for every case. The results from calculations 
v;ith the extreme temperature pairs thus defined the uncertainty in 
the RDMS M, P, and heat flux values. 


4.3 SODAR DATA ANALYSIS 


In Section 2.2 it was shown that under free convection conditions, 

2 

is related to the heat flux in this fashion: 


c2 3 , Va/3 -4/3 

S 4 <kg^ ^pCp) ^ 


(4.1) 


The result of taking the logarithm of this equation is a linear 

2 

relationship between log and log Z. 


log = log [| z 


(4.2) 


Note that the slope of this line is -4/3 and the intercept A is 


A 


log [| ( 


Xn2/3 

kg^ 


(_ 2 _)^/ 3 ] 
PCp J 


(4.3) 


Free convection is usually confined to a layer of the atmosphere 

several hundred meters deep. Thus, only a section of a log-log plotted 
2 

profile will have the characteristic -4/3 slope present if free 
convection exists (see Figure 4.9). A straight line with a -4/3 
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slope can then be fit Co that portion of the profile, and the value 

of the intercept is used Co produce a heat flux estimate. 

2 

Using the sodar data collected from every case, values were 

calculated and Chen averaged over each ten minute period from 900 to 

2 

1700 EST. From Chose averages, the log=log plots of versus height 

were produced. When a portion o' a profile showed a =4/3 slope, a 

'nomogram' of lines with a -4/3 slope and various intercepts calibrated 

2 

in values of heat flux was placed over the profile. Figure 4.10 

“2 

is an example of a 'nomogram' that goes from values of 5 to 60 W m 

-2 

by increments of 5 W m . The average heat flux for Chat ten-minute 
period based on Che free convection section of Che profile was 
calculated using the 'nomogram'. 
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5.0 RESULTS 

5.1 INTRODUCTION 

The results from this research need to be discussed in terms 
of tv;o separate measurement periods: 

1. The first three data sets collected in mid- and late July. 
(Hereafter referred to as the J cases.) 

2. The last four data sets taken during mid-September and 
early October. (Hereafter referred to as the S cases.) 

An examination of the heat flux estimates from both the sodar and 
the SDMS suggests significant differences bett'/een the tv;o periods. 

In the S cases, no usable sodar values of heat flux vjere obtained 
due to the absence of a clear-cut free convection regime, whereas 
the SDMS estimates of heat flux Increased beti^een the J and S cases. 
In the discussion beloii;, three primary causes for variations in 
surface heat flux between the J and S cases are evaluated: 

1. The change of ason. 

2. The aging of the crop between the two periods. 

3. The change in the soil moisture between July and September. 

5.2 RADIOMETRIC/MODEL DERIVED HEAT FLUXES 

Figure 5.1 displays the surface temperature ranges for every 
case. Between the J and S cases there is a drop of 6C and 7C 
in the average day and night temperatures, respectively. This step 
change is caused by the seasonal decrease in the amount and intensity 


Day I I - Night 



Figure 5.1 Day and night surface temperature ranges calculated for Cases 
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of oolar radlaclon. Solar radiation rocordo from the Ponn Stato 

°2 

Meteorological Oboervatory ohovjed 482 W m received on July 25 

=2 

(Case 3) verouo 422 W m on September 11 (Caoo 4), The scnith angle 

of the sun woo larger in the S caaoD, rooultlng in a docroooe in the 

intensity of the solar radiation. Although the temperature ranges 

for the S cases were cooler, the ROMS heat fluses (see Figure 5.2) 

= 2 

Increased an average of 25 W m beti/eon the J and S cases. Apparently, 
although the solar radiation was weaker in the S cases. It was 
partitioned in such a manner as to produce heat fluKes greater than 
those in the J cases. There is evidence that a reduction in the 
evaporative flux is responsible for this result. 

Along with the amount of rainfall, the age of the plants is an 
important factor determining the evaporation from a canopy. Vihen 
vegetation becomes senecent, reaching the end of its life cycle, 
evapotranspiratlon is greatly reduced. Observations indicated that 
senecent vegetation and a lack of rainfall were present in the S 
cases . 

Slroe the S cases occurred in the autumn, the crops at Rock 
Springs vi;ere brown and dying. To quantify the lack of rain in the S 
cases, an Antecedent Precipitation Index (API) was calculated for 
each case from the rainfall data at the Penn State Meteorological 
Observatory. The API, follovjing Blanchard et al., (1980)^, expresses 
moisture depletion as an exponential decaying function of the rainfall 
in this form. 


^Blanchard, B. J., McFarland, M. J., Schmugge, T. J., and Rhoades, E. : 
"Estimation of soil moisture with API algorithms and microxvave 
emissions." Unpublished manuscript. 
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API^ ° Pi 


(5.1) 


whoro 

API ° Gioioturo indoK 

P ° daily offoctive rainfall 

k ° doploeion conotant < 1 and a function of elmQ 
i ° day of tho ootlraato. 

An avoragG depletion conotant (k) for tho period of .920 uao uood and 

tho effective rainfall P v^ao related to the actual rainfall P„ in 

a 

L'.k8 faohion: 

P ° P (5.2) 

a 

The API values for each case are displayed in Figure 5.3. The average 
API in the S cases is one=half the value in the J cases. 

This change in soil and crop moisture was reflected in the ROMS 
measurements. Figure 3.4 is a block diagram of the ROMS H and P 
values for every case. Note that the average ROMS estimate of M 
for the S cases decreases by three tenths from the J cases. l-Jhen the 
moisture availability of a surface decreases, a larger amount of the 
net radiation is partitioned into sensible heat flux instead of 
evaporative flux. By examining the RDMS*-computed Bov^en ratios for 
each case (Figure 5.5), the change In the partitioning of the net 
radiation bett'feen the J and S cases is clearly demonstrated. 


J Cases 


S Cases 



Case Number 

Figure 5.3 The Antecedent Precipitation Index (API) 

for each case. An average depletion constant 
(k) of .92 was used. 









Figure 5.5 Bowen ratios calculated by 
the ROMS for Cases 1-7. 
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5.3 SODAR DERIVED HEAT FXUXES 

2 

When the sodar data were first analyzed, the profiles produced 

were 30 minute averages. Very few -4/3 slopes that extended over 100 

meters were observed In any of the cases. Figure 5.6 illustrates the 

state of the atmosphere when free convection exists. A matching layer 

where both surface and mixed layer scaling are valid Is present from 

a height of the absolute value of the Monln Length (L) to one tenth 

of the inversion height (.1 Z^) . A matching layer is most often 

found when winds are light, and solar heating Is strong. Thus, the 

2 

absence of extensive -4/3 slopes in the profiles was puzzling 
because of the clear skies and light winds (3-5 m/s) present when 
the data were collected. 

The original Intent of the 30-mlnute averages was to eliminate 

the effects of individual thermals and thereby show the average free 

convection layer. It appeared from the 30-mlnute averages however, 

that only a small portion, if any, of most 30-mlnute periods had free 

convection conditions existing. To see if free convection periods 

2 

were present on a shorter time scale, profiles averaged over 10- 
mlnutes were produced for all the cases. The 10-mlnute averages 
indicated that the sodar heat fluxes were also divided into J and 
S cases. With the sodar, however, the difference was much more 
drastic since the 10-minute averages improved the heat flux estimates 
in the J cases but failed to produce better results in the S cases. 

Although the J cases showed deeper free convection with the 10- 
minute averages, free convection was still present only intermittently 
throughout the day. The inconsistent presence of free convection in 


Free Atmosphere 


Z 


Mixed Layer 



Surface Layer 

77777 ////////////////////// 


Figure 5.6 


A schematic diagram of the atmosphere when free 
convection conditions exist. 
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che J caees is illustrafed by e»amlnaclon of two of the log=log plots 
2 

of versus height used to estimate the heat flux. Figure 5,7a Is 

2 

the average profile for 1350" 1400 EST on July 24 (Case 2), and 

Figure 5.7b is the average profile for the follov/lng ten minutes. It 

should be noted that the "4/3 slope line dravm on each plot is not 

the best linear fit to the -4/3 slope over the entire depth of the 

2 

free convection layer, but merely a line dravm through the value 

at 150 m to illustrate the -4/3 slope. Figure 5.7b has a free 

convection layer extending from 80 to 200 meters, while that for the 
previous ten-minute average indicates a shallov; layer of about 30 m 
at a height of 150 m. This lack of continuous free convection Is 
present throughout the J cases. 

In the S cases, the situation was even worse. Occasional free 

convection periods were noted in the early morning hours (900-1100 EST), 

2 

but were lacking in the remainder of the day. Figure 5.8 is a 

profile for 1400-1410 EST on September 11 (Case 4) showing the weak 
2 

C^ values characteristic of an S case afternoon. Thus, any heat 
fluxes that were determined in the S cases were weak and failed to 
display a diurnal trend. The lack of the expected diurnal variation 
increased doubts about obtaining any useful heat fluxes from the S 
cases' sodar data. 

Obviously something was disrupting the establishment of continuous 
free convection conditions at the site of sodar antenna. A mixed 
layer depth experiment using aircraft and sodar measurements had been 
performed at Rock Springs in 1978 (Lipschutz, 1978). Underwood (1981b) 
noted the presence of a low level inversion in the sodar data below 
the inversion height measured from the aircraft. In order to see if 





71 


a low level inversion was present in this OKperiment, the oiK hourly 
NEROS sounding data taken at Penn State during Cases 2 and 3 (July 24 
and 25) were examined. Figure 5.9 includes the series of soundings 
for Case 3 (July 25). Note the weak inversion at 920 mb in the 
afternoon (13Z) sounding that is still present in the evening (OOZ) 
sounding. Thus, even during the strong solar radiation of the J cases 
a lov; level inversion v;as present. 

Speculation on the cause of this inversion centered on the non- 
homogeneity of the terrain surrounding the site. As Figure 5.10 shows, 
the sodar is located on the east side of a broad valley. To the west 
are a ridge and elevated plateau. When the western ridge and plateau 
were heated, a plume of warm air at the surface of the plateau could 

have been advected by northerly winds at levels above the valley 

2 

floor . This elevated plume was a possible cause of the low level 
inversion plaguing the sodar results. 

The intermittent nature of the free convection conditions in the 
J case sodar data indicated that the low level inversion might 
periodically be disrupted by the penetration of a large thermal. 

At that time, short-lived free convection would exist and the sodar 

could estimate the heat flux. Figure 5.11 shows the quantitative 

false color display of the sodar data from Case 2 (July 24) and 

Case 6 (September 19). As shown by the Case 2 sodar returns, strong 

thermals were rising from the surface to about 200 m, whereas in 

Case 6 there was only weak thermal activity. The conclusion is 

that because of the reduced solar radiation in the S cases, the thermals 


Observed winds were northerly during every case. 
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Figure 5.11a Quantitative false color display of 
sodar data for Case 2. 
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Figure 5.11b Quantitative false color display of 
sodar data for Case 6. 
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were weaker Chan In Che J cases and were never able Co break Che 
inversion. Since a free convecCion layer couxd noc be escablished, 
no heac flujt esciraaces could be made from Che S cases* sodar daca, 

5.4 COMPARISON OF THE SDMS AND SODAR HEAT FLUXES 

Because Che sodar heac fluK values from Cases 4-7 v;ere considered 
unreliable, Che comparison beCween Che ROMS and sodar values was 
aCCempCed only for Cases 1-3. Figures 5.12-5.14 display Che RDMS heaC 
flux and solar radiacion escimaCes compared wlch sodar heaC flux and 
pyranomecer daCa, respecCively. 

The uncercalnCy in Che RDMS heac flux esCimaCes v;as caused by Che 
problems in deCerralning an area-averaged cemperature for Che heCero- 
geneous Rock Springs siCe. There were chree day/nighc CemperaCure 
pairs for each case ChaC defined Che esclmaced area-averaged value 
and ics error. The chree heac flux curves produced from Chose 
CemperaCui therefore defined the RDMS heat flux estimate and its 
uncertainty. In Figures 5.12-5.14, the ROMS heat flux is plotted as 
a swath with boundaries that lie halfway between the RDMS heat flux 
curve and the curves defining its uncertainty. The high and low 
heat flux curves themselves were not plotted as the error limits 
because those curves were often based on extreme values of M. For 
example, the low heat flux curve of Case 1 is calculated using an 
M value slightly greater chan 1. Thus, the most reasonable uncertainty 
for the RDMS heat flux was felt to be halfway between the middle 
curve and each extreme curve. 

The sodar heat flux values plotted in Figures 5.12-5.14 were 


corrected for three known error sources: 



1100 


1200 1300 1400 
Time (EST) 


1500 h 


L3 Heat flux and solar radiat 
and measured for Case 2. 
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Figure 5.14 Heat flux and solar radiation data calculated 
and measured for Case 3. 


82 


1. Uncercalncles In the calibration of the sodar. 

2. Excess attenuation from turbulence of the sodar signal. 

3. An Increase in the sodar signal caused by humidity fluctuations. 

The calibration process as described by Undenwod (1978) requires 
a measurement of the acoustic pressure at the antenna beam axis and 
the variation of output acoustic powet with RMS voltage given to the 
transducer. Both these measurements are determined in part using a 
microphone placed in the far field of the antenna beam pattern. If 
the microphone is not in the far field of the antenna, the sodar 
will be incorrectly calibrated. The far field of an antenna begins 
at about ten rimes the antenna diameter. Since the Aerovironment 
antenna is 1.2m in diameter, tie microphone should have been at least 
12m away. The calibration for this experiment was performed with the 
microphone approximately 7m from the antenna. Thomson (1981) estimated 
that the sodar heat fluxes could be in error from this problem by 5 %. 

Excess attenuation of the sodar signal is caused by turbulent 
spreading of the acoustic beam (Underwood, 1981a). The amount of 
excess attenuation is dependent on the antenna beam pattern and 
several atmospheric variables. Figure 5.15 developed by Underwood 
('1981c) was used to determine the excess attenuation term E in 
equation 2.1 for this experiment and was drawn using a friction velocity 
of 0.20 ms a frequency of 1600 hz, an antenna diameter of 1.2m, and 
an Inversion height of 1000m. By determining the average z^ ('v lOOOra) , 
z/z^ (.15), and z^/-L (200) for the first three cases from the RDMS 
results, E was estimated from Figure 5.15 as 10%. Since this value is 
only an estimate of the average E and the heat flux varies as the .75 
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power of C^, the reduction In the sodar heat flu» values due to excess 
attenuation was established at 10%. 

Humidity fluctuations can cause the sodar to overestimate the 
heat flux by Increasing the amount of backscattered energy. Coulter 
and Wesely (1980) derived a correction factor Yg^ that utilizes Bowen 
ratio measuiw ents to correct sodar heat flux values for humidity 
flu'-tuations. To a first approximation, for all values of the Boii^en 
ratio 


Ygg ~ (1 + O.O7/0g)~^ (5.5) 

where 6^ is the surface Bowen ratio. Multiplying Yg^ times the 

measured heat flux produces the correct value. Since no measurements 

of the evaporative flux were made in this experiment, an average RDMS 

value of the Bowen ratio in the first three cases was computed. Using 

this value for 6^ in equation 5.5, Yg^ was approximately 0.7. 

lue three aforementioned error sources were combined into a 

correction factor that lowered all the sodar-derlved heat fluxes by 15%. 

The corrected values are plotted in Figures 5.12-5.14. The error bar 

on each value represents a 10% uncertainty following Coulter and 

Wesely (1980). The actual comparison of the RDMS and sodar heat flux 

values was disappointing. In Case 3, the magnitudes were somewhat 

-2 

smaller, but there was a severe lack of agreement 30-50 W m ) in 
Cases 1 and 2. 

The disagreement of Case 1 (Figure 5.12) was characterized by the 

-2 

sodar values being about 45 W m lower than the RDMS heat fluxes 


'W" 
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bctv;een 1130 and 1330 EST. V/hlle the smooth ROMS heat flu» curve 

i-;as peaking, the sodar estimates suddenly dropped to a lower level. 

In the hours leading up to 1130 EST, the ROMS and sodar results 

-2 

disagreed in magnitude by 23 to 35 U m , but increased at the same 

_2 

rate. At 1330 EST, the sodar estimates gained 5 to 10 W m . Thereafter 
the rate of decrease of the sodar heat fluices matched that of the ROMS 
values. This afternoon period (1330-1700 EST) had better agreement than 
the morning period since the heat flux magnitudes differed by only 10 to 
30 W m"^. 

The sodar results of Case 2 (Figure 3.13) folloi-;ed the same basic 

pattern of Case 1. Between 900 and 1130 EST, the sodar values were 
-2 

35 to 45 W m too loii;, but increased at the same rate as the ROMS 

values. Between 1130 and 1500 EST, 4 out of the 5 sodar values v;ere 

-2 

at a level 10 to 25 W m below the values at 1105 and 1515 EST. The 

exception for this period was one sodar value at 1405 EST that was 
-2 

within 20 W m of the RDMS estimate. After 1500 EST, the sodar results 

-2 

were within 15 W m of the RDMS heat flux curve In magnitude and 
declined at the same rate. 

Case 3 (Figure 5.14) also had the same pattern of midday disagree- 
ment, but the RDMS and sodar va3.ues of the afternoon and morning 
periods were much closer than in the first two cases. The morning 

period of 900 to 1200 EST had two values within 10 W m ^ of the RDMS 

-2 

results. Three values were only 15 to 30 W m lower than the RDMS 
estimates and had the same rate of Increase. Between 1200 and 1400 KST, 

only two sodar estimates were obtained. These sodar values were 20 to 

-2 -2 
30 W m less than the RDMS values for that time and 10 to 15 W m 


below the sodar estimates at either end of the period. In the afternoon 

period for Case 3 (1415=1700 EST) , the ROMS and sodar magnitudes were 

-2 

close. All the sodar values v;ere within 10 W m of the ROMS values. 

Since an energy budget equation is the basis of the boundary 
layer model in the ROMS, if the model cannot accurately predict the 
amount of incoming solar radiation, the ROMS results will be invalid. 
Pyranometer data from the Penn State Weather Observatory were available 
for Cases 1-4 and 7. The model-predicted and measured solar radiation 
values are Included in Figures 5.12-5.14. The agreement betv^een the 
measured and predicted values is excellent. This result indicates 
that the ROMS heat flux curves are valid for Cases 1-3. The comparison 
for Case 7 is not as good, however (see Figure 5.16). The boundary 
layer model is overestimating the measured solar radiation by % 500 W m 
during the middle part of the day. In Case 7, therefore, the ROMS 
results may be incorrect since they are based on more solar radiation 
than was actually received. 

In the preceding discussion, a pattern of agreement in the trend 

of the values for a morning and afternoon period emerges. During 

-2 

the midday hours a period of strong (25 to 45 W m ) disagreement 
was present. Within this pattern, only Case 3 had comparable sodar 
and RDMS magnitudes for the heat flux. The boundary layer model and 
pyranometer solar radiation data agreed almost exactly for Cases 1-3. 



6.0 CONCLUSIONS AND RECOr-ffi-IENDATIONS 


I 


Tho purpooo of ehlo rosoarch v;ao to valldoto a method that uogd 
roraotely dorivod ourfaco tomporaturoo and a boundary layer model to 
ootimato Gurface enorsy budgoto and certain ourfaco parametoro. 

Bocauoe romotoly oonood temperature raeaourcraonto tend to roprooent 
valueo averaged over an area, the ourface flujjoo determined uolng 
Guch data muat aloo be area- aver agoo. To carry out on effective 
validation otudy. ground truth data muot be obtained from an Independent 
technique which aloo produces area=avoraged valuoo. 

One ouch method Is sodar. The oodar can produce heat flu:: 
estlraateo averaged over the same scale as the remotely oenoed tempera- 
ture data. By measuring the scattering of sound In a volume by 

temperature fluctuations using sodar, one can compute values of a 

2 2 

temperature structure constant, C^. C^ can In turn be related to 

the surface heat flu:c by a 4/3 power law when free convectlo** 

conditions e::ist In the atmosphere. 

In this research, the remotely sensed temperature data was 

obtained using a hand-held radiometer. A simultaneous collection 

of radiometric surface temperature and sodar data was made on seven 

occasions between mid-July and mid-October 1980. 

The comparison bettv-^en the radiometric data/model system (ROMS) 

and sodar derived heat fluxes proved disappointing. Free convection 

conditions were not present in Cases 4-7 and only intermittently In 

Cases 1-3. The sodar heat fluxes obtained in the first two cases 
-2 

were 25 to 45 W m lower than the ROMS values. In the third case, 

-2 

the sodar and RDMS heat flux values were ivichin 20 W m of each other. 
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The ROMS Id a relatively untootod method. It lo difficult to 
make strong concluoiono concerning Ito validity oinco the ROMS and 
oodar rooultD dioagroo. Hov/ever, the evidence preoentod Indlcateo 
that the oodar valueo could be in error by greater than the 102 
figure given by Coulter and V/eoely (1980). A pooolblo oourco for 
the poor oodar reoulto vmo the preocnce of a pcrolotont lov; level 
Inveroion. Thio Inverolon v/as formed by a warm plume advocted over 
the valley oodar oite from an elevated heated plateau. 

The results do ouggeot, however, that the true heat fluKOo v;oro 
not underestimated by the ROMS. Therefore, the Bovjen ratios over 
well-watered vegetation are likely to be quite small. Computation 
of an Antecedent Precipitation Index (API) suggests a decrease in 
the soil moisture between July and September. This decrease is also 
present in the RDMS results. 

Any future attempts to validate the remotely-sensed surface 

temperature and boundary layer model method should take heed of the 

difficulties encountered in this research. If the radiometric 

surface temperature data is to be used again, a less subjective 

method to determine an area-averaged temperature s.iould be developed. 

Steps should also be taken to improve the quality of the sodar results 

by insuring that the effects of complex terrain are minimized. In 

addition, the technique successfully used by Coulter and Wesely (1980) 

to obtain sodar-derlved heat fluxes should be emulated. Their 

procedure Includes simultaneous heat flux measurements using the 

2 

eddy correlation method and real time display of C^ profiles to 
determine conclusively If free convection exists. Improving the 


sodar resulte and uolng bettor surface temperature data will enable 
atronger conclusions to be dravm about the validity of the remotely 
sensed surface temperature and boundary layer model method. 
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APPENDIX A 

SURFACE TEMPERATURE ERROR FROM ASSUMING A GROUND EMISSIVITY OF ONE 


Many investigators (Fuchs and Tanner (1966), Conway and Van Bavel 
(1967), Fuchs and Tanner (1968)) using radiometers to determine 
vegetation and surface temperatures, estimated the surface emissivity 
and atmospheric dovmwelling radiation in order to increase the 
accuracy of their measurements. An investigation into the error 
caused by assuming a surface emissivity of one is shown below. From 
that explanation, it was determined that for the purposes of this 
research it was not necessary to take the reflected dotmt^elllng 
radiation into account. 

Figure A.l shows that the longwave radiation reaching the 
radiometer is a combination of ground and reflected downwelllng 
radiation from the atmosphere. This situation can be written 
mathematically using Planck's Law as 

aT^ = e aT^ + (1-e Je a T ^ (A.l) 

e go g a a 


where 

a = Stefan Soltzman constant, 

T^ = effective temperature as measured by the radiometer, 

e = emissivity of the ground, 

8 

T^ = actual ground temperature, 
e = emissivity of the air, 

T^ = temperature of the air. 
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Figure A.l Illustration of the longwave radiation 
reaching a radiometer, 

G = emissivity of the ground 
8 

° emissivity of the air 
temperature of the air 
= temperature of the surface 
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Solving aquation A. I for the actual temperature yields 


a 

O 


a T 




a 


e o 
g 


a T 


1/4 


(A.2) 


By using some reasonable values in (A.2), an estimate of the 
difference bett>jeen the real and actual values for this research can 
be obtained. It if is assumed th..<t 

the effective temperature T is 300 IC, 

6 

the air temperature T is 300 K, 

the emissivity of the air e^ is 0.8 (Sellers (1967)), 

and the emissivity of the ground e is 0.98 (Fuchs and 

g 

Tanner (1966)), 

then the actual ground temperature is 300.31 K. Thus, the actual 
temperature is only .3C warmer than the temperature measured by the 
radiometer. It is clear then, that to achieve an accuracy of ± O.IC, 
one must determine the actual emissivity of the ground. 

However, two other error sources in this experiment overwhelm 
the emissivity error. 

1. Tanner et al. (1967) noted a ± 0.3C variation in time 

with their temperature measurements just from air turbulence. 

2. A pre-experiment calibration of the radiometer showed that 
it consistently overestimated the temperature of a stirred 
water surface by 1.5 to 2C. This error was accommodated by 
subtracting a constant 1.7C from the measurements. 


The additional effort of determining the surface emissivity of 
the heterogeneous cropland at the Rock Springs site is clearly not 
worth the effort. The error from assuming a surface omissivlty of 
one is lost in the radiometer and turbulence errors which combined, 
make each measurement accurate to ± 0.3C. Since satellite radiometric 
and the boundary layer model temperatures are accurate only to £ 1 or 
2C, the error of t 0.3C is acceptable. 
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APPENDIX B 

SURFACE TEMPERATURE MEASUREMENTS FROM AN AIRBORNE RADIOMETER 


In Case 6, In addlclon to the regular surface temperature 
meaeurements , temperature data v;ere collected from an airborne 
radiometer. An Everest Infrared Thermometer v/lth a 3 degree field 
of view was flovm four times at about 150m over the Rock Springe 
site. During the overflights, temperature values for various crops 
were recorded. The raw data are presented in Table B.l. The Everest 
radiometer does not assume an emlsslvlty of one in producing its 
temperature values. Unlike the Barnes radiometer used on the ground, 
the Everest user specifies an emlsslvlty from 0 to .99. Thus, the 
theory discussed in Appendix A must be altered to Include the radiometer 
radiometer emisslvity e^. Putting in equation A.l we obtain 


a T 


a T + e a T (1-e ) 
o a a g 


(B.l) 


We want to calculate how much the radiometer-measured (effective) 

temperature T changes v;hen e is changed from a value of one. If 
6 r 

we assume e is constant and known and differentiate (B.l) with respect 

5 

to Gj, we get 


3 3T . 

S IT ^ ^ ° 


3T 


Solving (B.2) for ■5— , the result is 

O 


(B.2) 


Pass 
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(B.3) 


Uolng finite difforGncGo,, we can qoIvq for the change in T cauoed by 

G 

a change in from a value of one. 

T 

AT ° (-r) Ae (B.4) 

e 4 r 


To compare the Everest temperature meQOuremento with the tempera^ 
tures made on the ground* a correction in the Evereot T^ for an of 
.93 must be made following the theory prevlouoly dlecuooed. If the 
Everest measures a temperature of 300 K using an of ,98, Ae^ io 
-0.02 and the corrected Everest temperature is 1.3 IC warmer than the 
measured value. 

The average Everest temperature for the site was 24. 3C. Adding 
in the correction factor of 1.5C, the result is 23. 8C. This value is 
1.3C cooler than the 27 3C estimated from the ground data. However, 
the uncertainty of the ground estimate is ± 2.3C and no water vapor 
correction was made to the Everest value. The biggest difference 
betx^een the two sets of temperatures was the spread in the values. 

The ground data v;ent from 21. 1C to 29. 6C v;hile the Everest values 
ranged from 23. 6C to 28. 3C. The Everest values narrow dot-m, but do 
not extend beyond the range of ground values. The agreement between 
the estimated area-averaged values Indicates that the subjective 
analysis performed on the ground temperature data produces a reasonable 
area- averaged temperature value. 
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APPENDIX C 

DETERMINATION OF RDMS HEAT FLUX FOR CASE 3 

Tho boundary layer modol, after having been ouitably Initiallsod 
for Caoo 3 (July 23, 1980), wao run for the period 0600 to 2400 EET 
16 tloiQO uaing different moioturo availability (M) and thoriBal 
inertia (P) combinationo that spanned thooe pararootero likely rangoo 
(M o 0.23=1, P ° 0,003=0.08), The ourface temperoture data for 
Caoo 3 v;ao collected at 1300 EST and 2400 CST. At thooe tlineo, the 
model output for each M and P combination mo extracted and stored. 

AO a olmple t-^ay of repreoenting the 16 oeto of model output, oocond 
order regression equations were formed vjlth surface temperatures ao 
predictors and M, P, or a face flu» as the predictand. The 
regression equations had the form 

X = Cq + CjTj + . Cj/ 

v;here X is M, P, heat flux, etc., are the regression coefficients 

for the parameter, T^ is the day surface temperature, and T^ is the 
night surface temperature. In this research, the Important equations 
vtere those predicting M and P and the regression coefficients for 
those equations are given in Table C.l. 

The next step was to use the observed surface temperatures to 
predict an M and P value for this case. T^ and Tj^, estimated at 
30. 3C and 14. 3C, respectively, produced an M of 1.0 and a P of 0.068. 

To determine the RDMS heat flux estimate for Case 3, the model was 
run again v/lth the same set of initial conditions, but M and P v;ere 
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hold fiisod at liO and 0.068, roopoctlvoly. Tho rooulting diurnal 
hoae flujs curve and 1 to uncoreainey aro plotCod ao a owath In Figure 

5.U. 


Table C.l Rogroooion coofficionto for M ond P in Cooo 3. 



M 

P 

s 

4.75 

0 *446 

^1 

“0.165 

=0.0201 


0.0376 

=0.00504 


0.00147 

0.00231 

Ca 

=0.00289 

0.000451 


